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Abstract 

We investigate baryogenesis in the i/MSM, whicli is tlie Minimal Standard Model 
(MSM) extended by three right-handed neutrinos with Majorana masses smaller than the 
weak scale. In this model the baryon asymmetry of the universe (BAU) is generated via 
flavour oscillation between right-handed neutrinos. We consider the case when BAU is 
solely originated from the CP violation in the mixing matrix of active neutrinos. We 
perform analytical and numerical estimations of the yield of BAU, and show how BAU 
depends on mixing angles and CP violating phases. It is found that the asymmetry 
in the inverted hierarchy for neutrino masses receives a suppression factor of about 4% 
comparing with the normal hierarchy case. It is, however, pointed out that, when ^13 = 
and 023 = baryogenesis in the normal hierarchy becomes ineffective, and hence the 
inverted hierarchy case becomes significant to account for the present BAU. 



1 Introduction 



The origin of the baryon asymmetry of the universe (BAU) is one of the most mysterious 
problems in particle physics and cosmology, since the Minimal Standard Model (MSM) and the 
Big Bang cosmology cannot answer it. So far various mechanisms for generating BAU have 
been proposed [1]. One promising possibihty is the so-called leptogenesis scenario [2] (see also 
Ref. [3]), where the origins of neutrino masses as well as BAU can be addressed at the same 
time by introducing right-handed neutrinos with superheavy Majorana masses. The observed 
tiny masses of neutrinos can be naturally understood by the seesaw mechanism [1]. Further, 
the lepton asymmetry generated by decays of right-handed neutrinos can be a source of BAU. 
In the simplest thermal leptogenesis, the required Majorana masses is heavier than about 10^ 
GeV [5]. 

It should be, however, noted that the connection between the origins of neutrino masses and 
BAU can be obtained even when Majorana masses are below the weak scale. One interesting 
possibility is the so-called z/MSM [6l [7], which is the MSM extended by three right-handed 
neutrinos with masses smaller than the weak scale. In this model the problems of neutrino 
masses, BAU and also dark matter can be solved simultaneously. One attractive advantage 
of the i/MSM lies in the fact that the direct detection of right-handed neutrinos becomes 
possible [S]. 

In the i/MSM BAU can be generated by invoking the mechanism via flavour oscillation of 
right-handed neutrinos [9]. (See also Ref. [HI [10].) In this mechanism the lepton asymmetry 
is not generated for temperatures of interest because of the smallness of Majorana masses, 
which is very different from the leptogenesis scenario. The lepton asymmetry is separated into 
left-handed and right-handed leptonic sectors due to the CP violations in the production and 
oscillation of right-handed neutrinos. Then, the asymmetry stored in the left-handed sector is 
partially transferred into the baryon asymmetry through the rapid sphaleron transitions 

One of right-handed neutrinos in the i/MSM, which is a candidate of dark matter, plays no 
essential role to generate BAU as well as neutrino masses observed in the oscillation experiments, 
since its Yukawa coupling constants should be very suppressed. The rest two are responsible 
to BAU via their flavour oscillation, but also induce the masses of active neutrinos through 
the seesaw mechanism. Therefore, physics of these two right-handed neutrinos connects BAU 
with the neutrino parameters of active neutrinos, i.e., mass hierarchy, mixing angles, and CP 
violating phases. 

In this letter we would like to extend the analysis in Ref. [B]. Under the considering situ- 
ation there are three CP violating phases in the leptonic sector which can be a source of the 
asymmetry. Especially, we concentrate here on the case when BAU is originated only from 
the CP violation in the mixing matrix U of active neutrinos, namely, the Dirac phase 6 and 
Majorana phase rj in U. We then present the analytical expression of BAU shows explicitly how 
BAU depends on these CP phases and the mixing angles of active neutrinos. Moreover, we also 
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perform the numerical estimation of BAU and justify the vahdity of the analytical expression. 

2 The z/MSM 

We begin with the brief review of the model under consideration, i.e., the z/MSM [SI E]- It 
is the MSM extended by three right-handed neutrinos uju (J = 1,2,3), which Lagrangian is 
given by 

^uMSM = ^MSM + iWll^ Vri - [Pal $ ZZ/j/ + ^ Vjfj Vrj + /l.C.j , (1) 

where £msm is the MSM Lagrangian. F^i are neutrino Yukawa couplings, and $ and 
[a = e, fi, t) are Higgs and lepton weak-doublets, respectively. The Majorana masses of right- 
handed neutrinos are denoted by Mj which are taken to be real and positive without loss of 
generality. Here and hereafter we work in a basis in which the mass matrix of charged leptons 
is diagonal. In this model, neutrinos also obtain the Dirac masses, [Md]^/ = Fq,/($) (($) is a 
vacuum expectation value of the Higgs field), after the electroweak (EW) symmetry breaking. 

The distinctive feature of the model is the region of the parameter space of Eq. ([1]), i.e., we 
restrict ourselves in the region 

|[M,,],,|<M, <100GeV. (2) 

In this case the seesaw mechanism works, and mass eigenstates of neutrinos are then divided 
into two groups. One group consists of active neutrinos z/j (i = 1, 2, 3). Their masses are found 
from the seesaw mass matrix = —M^Mf^M]^ as 

M^U* = dicig{'mi,m2,ms) , (3) 

where U is the mixing matrix of active neutrinos. The other one consists of sterile neutrinos 
Nj which are almost the right-handed states Nj ~ upu, and their masses are approximately 
given by the Majorana masses Mj. We then find the neutrino mixing as 

i^La = Uai Vi + e„j N\ , (4) 

where 9^/ = [M]:,]ai / Mi are the active-sterile (left-right) mixing matrix. We should stress 
that sterile neutrinos Nj here are originated from the right-handed neutrinos in the seesaw 
mechanism. Thus, we simply say Nj as right-handed neutrinos from now on. 

In the z/MSM three right-handed neutrinos play important roles in cosmology. One of them, 
say A?"!, is a candidate for dark matter of the universe |7j. It is beyond the scope of the present 
work to describe this issue. However, one thing being important for later discussions is that 
the Yukawa couplings of A^i should be highly suppressed to realise a successful dark matter 
scenario. (See the details, e.g., Refs. |12l [TOl [131 ^ result, the contribution from Ni 

to the seesaw matrix My becomes negligible [7]. Furthermore, as shown in Ref. [6], A^i plays 
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essentially no role to generate BAU. Therefore, we take F^i = for simplicity in the rest of 
this analysis. 

The other right-handed neutrinos, N2 and N^, are then responsible to the masses and 
mixing angles of active neutrinos. Notice that the lightest active neutrino becomes massless in 
our approximation. Further, the flavour oscillation between and N3 in the early universe 
can be a source of BAU through the mechanism proposed in Ref. [9], as we will show below. 
In the i/MSM, therefore, BAU is related to the parameters of active neutrinos through physics 
of A^2 and N3. 

The neutrino Yukawa matrix F for N2 and N^, which is a 3 x 2 matrix, can be expressed 
without loss of generality as [I5] 



UDl'^QDf. (5) 



Here parameters of active neutrinos are their masses Dy = diag(mi, m2, ^3) and the mixing 
matrix 



C12C13 S12C13 si3e 



U = -C23S12 - S23Ci2Si3e'^ C23C12 - S23Si2Si3e*'' S23C13 X diag(l , e*'' , 1) , (6) 



S23S12 - C23Ci2Si3e*'' -S23C12 - C23Si2Si3e*'^ C23C13 

with Sij = sin 6ij and Cij = cos 6ij. The Dirac and Majorana phases are denoted by 6 and rj, 
respectively. Since we set F^i = 0, the masses of active neutrinos are 

^3 = fnatm > ^2 = ""^soi > =0 in the NH case , 

^2 = \j "^atm + "^soi > ^1 = ^/^^^ > nis = m the IH case , (7) 

The observational data of mixing angles are s^2 = 0-318^o;o48, s^^ = O.SOIq'iL ^13 — 0.053, 
respectively, and masses are m^^j = Am^i = {7.59^q^q) x 10~^ eV^, and m^tm = I'^^^ail — 
(2.40+ 

ass) X 10 eV^ (at the 3a level) [I6]. Hereafter, we shall adopt the central values unless 
otherwise stated. 

On the other hand, parameters of N2 and are their masses -Dat = diag(M2, M3) and the 
3x2 matrix 



= I cos tu — sin tu I in the NH case 
^ sin w ^ cos u 



cos u — sm a; 

fi = I .^sina; .^cosa; | in the IH case , (8) 


where ^ = ±1 and u is an arbitrary complex number. 

In the considering situation there are three CP violating parameters, 6, 1] and Imu, in 
leptonic sector, which can potentially contribute to the generation of BAU. In this analysis we 
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concentrate on the case in which BAU is originated solely from the CP phases in the mixing 
matrix of active neutrinos, and find the dependence on 6 and rj as well as the mixing angles 6ij 
by taking Imoj = 0, 
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3 Baryogenesis via Neutrino Oscillation 

Let us then discuss how BAU is generated in the i/MSM through baryogenesis via oscillation 
of right-handed neutrinos [9]. In the considering model the lepton-number violations due to 
Majorana masses are ineffective for high temperatures T > 100 GeV [H]. In order to generate 
the baryon asymmetry, thus, it is crucial the lepton asymmetry is distributed into left-handed 
leptons La and right-handed neutrinos Nj for the temperatures of interest rather than its 
creation. Then, the asymmetry stored in the left-handed sector is partially transferred into 
the baryon asymmetry due to the B + L breaking sphaleron transition which is rapid for 
T>Tw^ 100 GeV [II]. 

We denote the asymmetries of numbers of A'"2,3 and by AA^2,3 and AL^, respectively. 
These asymmetries are estimated by solving the kinetic equations for their density matrices 
Pnn and PlJ*§, which are given by 

= [Hnn + Vn, Pnn] - ^{r^v, PNN - Pnn) + ^-^y^TF\pll - Pll)^ ' (9) 
= [HlL + VL,PLL]-'-{rL,PLL-pl',}+'^TF{p^^-p';^^)F^, (10) 

where note again that F is the 3x2 matrix for neutrino Yukawa couplings of N2 and N^. Hp^j^ 
and H^j^ denote the Hamiltonian when F = 0. The effective potentials and the destruction 
rates for and La are 

Vn = ItF^F, rr, = 2s^Vr,, (11) 

o 

Vl = ^TF^F, rL = 2s^VL. (12) 
lb 

where s,^ ~ 2 x 10"^ [9]. Notice that these expressions are valid as long as T is sufficiently 
higher than T^r. The kinetic equations of the density matrices for the anti-particles p^^y ^ind 
Pll are obtained by replacing F — )• F* in Eqs. and (ITUl) . 

These equations include the medium effects of surrounding hot plasma, i.e., the thermal po- 
tentials which describe the coherent oscillations of right-handed neutrinos induced by Vjy, and 
the decoherent terms which describe the production and destruction of Nj [9]. Furthermore, we 
include the terms which express the exchange of asymmetries between left and right-handed sec- 
tors [6]. It is found from Eqs. (E]) and (JIO]) that AA'tot + ALtot = 0, where AA'tot = E/=2,3 



The impact of Imo; on BAU will be discussed in elsewhere |17| . 
"^"^ The density matrix denotes the sum of Pulvl Pclblj which are the same in the temperatures under 
consideration due to SU(2) l symmetry. 
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and ALtot = Sa=e/xr^-^a' which is crucial in the considering baryogenesis scenario as men- 
tioned above. For further details of these issues, please see Ref. [6]. 

The coupled equations ([9]) and (fTOj) can be solved not only numerically, but also analytically 
by using the perturbative expansion of the Yukawa coupling constants F |6]. The initial con- 
ditions are taken as Pnn{S^) = Pnn{^) = and Pll{0) = PLli^) = Pll- Then, we can estimate 
the asymmetries as ANj = [pnn]ii - [Pnn]ii and ALa = [pLL\aa - [pillaa- In the following 
we will present the analytical expression of the active flavour asymmetry as well as BAU at the 
leading order of F . We will also show the results from the numerical solutions of Eqs. ^ and 
( 1T0|) . which confirm the validity of the analytical expressions of the asymmetries. 



4 Active Flavour Asymmetries 

First of all, we discuss the yield of the active flavour asymmetry ALa (a = e, p, r). The leading 
order contribution to ALa is induced at 0{F^) and the analytic expression at the temperature 
T can be written as [6] 

ALa{T) = ^A'^,^MTjT), (13) 

L 

where Mq ~ 7.1 x 10^^ GeV and the CP asymmetry parameter is defined by 



^32 = Im 



The typical temperature of the generation of ALa, T^, is given by 



where AM^2 = — M|. The evolution of ALa is described by the function J32 

■^32(3;)= / dxi / dx2sm{xl — xl) . (16) 
Jo Jo 

It is found that J32(x) ~ ^x^ for a; <^ 1 while it is oscillating around and approaching to the 
value 

21/3^3/2 

J32(a:)U»i = ^^;^ = 0.69. (17) 

Therefore, the active flavour asymmetry scales as ALa oc for T ^ Tl and takes a constant 
value for T <^ as 

181(5/6) "^^2 (AM|2)2/3 • ^^^^ 
We can see that ALa for T ^ T;, is enhanced when N2 and iVs are quasi degenerate [6] . 
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Now we would like to express the CP asymmetry parameter by using the parameters 
of the neutrino Yukawa couplings in Eq. ([5]). In the NH case we evaluate 

.a c ■ OT) M3M2(m3m2)^/^(m3 - 7712) 
Ann = f sm ztieoo 1 x 



3/2 ^1/2 
4($)^ 



e sin 2Rec. ^^^^^l^t"^-' x S^^ x a„ , (19) 



where is 



Snii, — 1 ''"m ; (20) 

where Vm = rnsoi/rria.tni- The parameter Oq can be evaluated as 

Oe = — sin 6*12 sin 26^13 sm{6 + rf) , 

= + sin^ 6^23 sin 6^12 sin 26*13 sin (5 + r]) — sin 26^23 cos 6^12 cos Oi^ sin 77 , 
ar = + cos^ 6*23 sin 6*12 sin 26*13 sin(5 + r/) + sin 26*23 cos 6^12 cos 6^13 sin rj . (21) 

We find that the active fiavour asymmetry depends on parameters of right-handed neutrinos 
as ALq, oc ^sin2Rea;M2M3/(AM|2)^/^, and hence vanishes when Recu = 7111/2 [n is integer). 
This is simply because the flavour oscillation between N2 and N^, which is essential to the 
considering mechanism, disappears. On the other hand, the dependence on the mixing angles 
and CP violating phases of active neutrinos is summarised in a^. 

As already pointed out in Ref . [B] , the total asymmetry of active flavours vanishes ALiot = 
at the leading 0{F^) since ^32 = Im [F^F]23 [^"^^132 = 0, which can be seen as + + 

= in Eq. (12T|) . It is interesting to note that the active flavour asymmetries depend on CP 
violating phases in two ways, i.e., sin(5 + rj) and sin 17, and ALg only depends on the former 
one. We also flnd that ALg vanishes and AL^ = — AL,- when 6*13 = 0. In this case = and 

= —Or = — sin 2^23 cos 6*12 sin ?7, and the asymmetries only depend on Majorana phase r] as 
expected. 

On the other hand, in the inverted hierarchy case, the CP asymmetry parameter ^432 can 
be written as 

.a t • 00 M^M2{m2miY/'^{m2-mi) 
Ann = f sm 2Rea; ^ —r x Oq, 

= e sin 2Rec. ^^^^"^'^t"^'" - Sm. x a„ , (22) 



4($) 
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where is found as 



+ sin 26i2 cos^ 6^13 sin 77 , 
- sin^ ^23 cos^ 6^13 + cos 2^23 



+ sin 26'23 sin 6 



'13 



sin 2^12 sin 
cos rjsm.5 — cos 26'i2 sin 77 cos 5 



COS^ 6*23 COS^ 6'i3 



cos 2^^ 



23 



sin 2^12 sin rj 



— sin 26'23 sin 
Further, we have introduced 



cos 7] sin 5 — cos 2^^i2 sin r] cos 5 



1 - ^ + 0{rl) 



(23) 



(24) 



We find that the active fiavour asymmetry in the IH depends on the parameters of right- 
handed neutrinos in the same way as in the NH. The total asymmetry of active fiavours at 
0{F'^) vanishes ALtot = as in the NH case. 

It should be noted that, comparing with the NH case f|T9l) . the IH case receives of 

suppression factor of Sm^ ~ 0.04, and thus the production of active fiavour asymmetries in IH is 
less effective than NH apart from a^. Further, ALq, depends on CP violating phases differently 
from the NH case. When ^13 = 0, we find that 



+ sin 2^^i2 sin rj , 

— sin 26^12 cos^ 6^23 sin r] , 

— sin 26*12 sin^ 6*23 sin 77 . 



(25) 



Thus, as in the NH case, the asymmetries only depend on Majorana phase 77, as expected. In this 
case, however, ALe does not vanish even when 6*13 = 0, and we have AL^ = AL^- = — 0.5ALe 
for ^23 = vr/4. 

The total asymmetry ALtot is induced at 0{F^) and ALtot = — AA^tot- It is found in Ref. [6] 
that 



ALtot (T)= - 



S<f> 



1 /=2,3 a=e,/i,T 



1=2,3 a=e,ii,T 



M3 



n 



where the evolution of ALtot is described by the function K: 



K{x) 



dxi J32ixi 



(26) 



(27) 



We find that K{x) ~ x^/40 for x ^ 1 while K{x) ~ 732(00) x for x ^ Icj Then, we normalise 



We numerically find that K{x) can be fitted as K{x) — ax ~ b ioi 2 < x < 100, where a — 0.68 ~ ^32(00) 
and 6 = 1.4. 
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Figure 1: Function k,{x) in Eq. (125]) . 



K{x) as 



K{x) = 732(00) 



(28) 



where the behaviour of the function k{x) is shown in Fig. [H 

It should be noted that we have so far assumed that N2 and are out of equilibrium in 
order to avoid the wash-out of the asymmetries. When the interaction rate of right-handed 
neutrinos (1TT|) is F^r < H (the Hubble parameter) till the sphaleron transitions are switched 
off at T = Tw — 100 GeV, the eigenvalues of F'^F should be smaller than 2.8 x 10~^^. This is 
translated into the upper bound on the mass of right-handed neutrinos as M2,3 ^17 GeV [6J 
for both the NH and IH cases by taking N2 and are quasi-degenerate, which is required 
to explain the observed BAU as we will show. We thus restrict ourselves in this mass region 
below. 

5 Baryon Asymmetry of the Universe 

Now we are at the position to present the analytical expression of BAU and study how it 
depends on masses, mixing angles and CP violating phases of active neutrinos. The baryon- 
to-entropy ratio, Yb = ub/s {ub and s are the baryon and entropy densities at the present 
universe, respectively), is obtained from ALtot as 



where Cs = 7.0 x 10"'^. Notice that only the asymmetry in the left-handed leptons is transferred 
via the sphaleron processes into the baryon asymmetry as AB = —|| ALtot [IS], and this 
transition is switched off for T < Tyy. 

We take the typical temperature f lTSj) for the generation of AL^ as > Tw, and A^2 and 
A^3 are quasi-degenerate, i.e. AM < Mn where = Mn + AM/2 and M2 = Mn - AM/2. 



Yb 



Cs — ALtot (TV) 5 



(29) 
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These are arranged in order to enhance the production of the baryon asymmetry. In this case, 
we find that the analytical expression of Yb for both the NH and IH cases is given by 



Yb = 1.8 X 10-" k{Tl/Tw) ° ^ "^atm"W ^ ^ 



4.7X10-.CP(— )(^) (^^^) . (30) 



where we have taken k{Tl/Tw) = 1 in the last equality. The CP asymmetry parameter 5cp for 
BAU is expressed as 

^cp = i sin 2Rea; ■ Sm, 5^ . (31) 

Here Sm^ is the parameter concerning with the neutrino mass hierarchy in Eq. fl20p or (1241) 
for the NH or IH case, respectively. Therefore, Yb in the IH case is suppressed by about 4% 
compared with the NH case apart from 5ij. 

It can be seen that to get the sizable Yb the mass degeneracy of N2 and N^^ at a rather 
high accuracy is required. Note, however, that such a small mass difference is stable against 
the radiative corrections due to the smallness of Yukawa coupling constants of neutrinos under 
consideration. We also find that a larger Yb can be obtained for a larger M^r as long as 
^ 17 GeV. Otherwise, N2 and get in thermal equilibrium and the asymmetries are 
washed out. It should be noted that Yb is proportional to ^sinRecj, which means that the 
sign of Yb cannot be uniquely predicted even when all the parameters of active neutrinos are 
experimentally determined. 

In this expression we have introduced the CP asymmetry parameter 5i, to describe how Yb 
depends on the mixing angles and CP violating phases in the mixing matrix U . The analytical 
expression for 5^ for both the NH and IH cases are presented in Appendix |Al From now on 
we shall study these dependence by using the analytical expression as well as the numerical 
estimation of Yg, which is obtained by solving numerically Eqs. ^ and flTU]) . 

In the NH case we find that the leading term of 5^ in the expansion of is given by 

^u = ^ sin ^12 sin26li3[cos^^i3(3 + cos 4^23) - 4sin^^i3] sin (5 + rf) 

+ cos 6*12 sin 46^23 cos^ 6^13 sin rj + C(rm) • (32) 

It is seen that 5^, depends on the CP violating phases in two ways, i.e., 5 + rj and 77. 

We then find that when the mixing angle in the atmospheric neutrino oscillation is maximal 
^23 = 7r/4, the leading term of 5y depends only on the sum of the CP phases 5 + rj. This 
behaviour can be understood in the left panel of Fig. [21 which represents the parameter region 
accounting for the present observation data of Yb in the 5-7] plane obtained by the numerical 
solutions of Eqs. ([9]) and (1101) . When 6'23 is slightly smaller than 7r/4, the allowed region becomes 
drastically changed as shown in the right panel of Fig. [2l Thus, the deviation of 6*23 from 7r/4, 
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8/ji 5/7t 



Figure 2: Parameter regions in the 6-7] plane accounting for the observed baryon asymmetry in 
the NH case are shown by hues. The red sohd hues are for ^ = +1 while the blue dashed lines 
are for ^ = -1. Here we take M3 = 5 GeV, AMiJMl = 10-^ Reu = n/A, and sin^ = 0.053. 
We take sin^ 6^23 = 0.5 (left) and 0.36 (right), respectively. 

which will be tested in future oscillation experiments, is significant to determine Yb in the 
considering case. 

On the other hand, when ^13 = 0, the CP asymmetry parameter becomes 



including the higher order term of r^. In this case the asymmetry depends only on the Majorana 
phase 77 as expected (since the Dirac phase S always appears together with S13). It is very 
important to note that 6^ = when ^13 = and 623 = it /A. In this case, the generation of BAU 
in the NH case is ineffective and Yb at 0{F^) vanishes. 

Next, we turn to consider the IH case, where the CP asymmetry parameter 6^, at 0{r^) is 
estimated as 



It is then found that 5i, at the leading order depends only on the Majorana phase when 623 = 
7r/4, which should be compared with the NH case. This behaviour is shown in the left panel of 
Fig. [31 Moreover, we discover that 6^ in the IH case does not vanish even when 623 = it /A and 



5u = sin 4023 cos ^12(1 - rm cos^ 612) sin 77 , 



(33) 




- sin 26^12 cos^ ^13 [ — 5 — 3 cos 4^23 + cos26'i3(7 + cos 4^23)] sin 
+ sin 46*23 cos^ 6*13 sin 6*13 (sin 5 cos — cos 26*12 cos 5 sin 77) + C(r^) . 



(34) 



0; 



6, = ^[1 + {1 + r^)^/2 + 3[1 - (1 + r^)^/2] cos 2^12] sin2^i2 sin 77 
= sin 2^12 sin 77 1 + ^(1 -3cos2^i2)r^ + C(r^) . 



(35) 
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Figure 3: Parameter regions in the S-rj plane accounting for the observed baryon asymmetry 
in the IH case are shown by shaded regions. The red sohd hues are for ^ = +1 while the blue 
dashed lines are for ^ = -1. Here we take M3 = 15 GeV, AM^^/Ml = 10"^, Reu = n/4, and 
sin^ ^13 = 0.053. We take sin^ 6^23 = 0.5 (left) and 0.36 (right), respectively. 

This is one important feature of generating BAU in the IH case. 

Therefore, in the IH case is smaller than the NH case in the most of the parameter space, 
since it receives extra suppression factor of in Eq. (^^. See Fig. HI However, when ^13 is 
very small and also 623 is almost maximal, baryogenesis in the NH case becomes ineffective and 
hence the IH of active neutrino masses is the essence of accounting for the present BAU in the 
uMSM. 

6 Conclusions 

We have studied baryogenesis in the z/MSM via flavour oscillation between right-handed neutri- 
nos N2 and A''3. In particular, the case when BAU comes solely from the CP violating phases in 
the mixing matrix of active neutrinos has been investigated. We have presented the analytical 
expressions of BAU for both the NH and IH cases of active neutrino masses, and have demon- 
strated how the present value of BAU depends on the Dirac and Majorana phases as well as 
the neutrino mixing angles. We have shown that BAU depends on the neutrino mass hierarchy 
and Yb in the IH case receives the suppression factor of Sm^ ~ 4%. It has been found that Yb 
is very sensitive to the mixing angles 623 and 6*13. When 623 = 7r/4, the leading contribution to 
Yb is proportional to sm{S + rj) for the NH case while to sinry for the IH case. Moreover, when 
623 is almost maximal and 613 is very small, the CP asymmetry parameter in Yb vanishes and 
no baryon asymmetry is generated (at least the leading 0{F^) contribution) in the NH case. 
In this case, the IH case is required to explain the observed BAU. 
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Figure 4: in terms of sin6'i3 are shown by lines. The red sohd hne is for the NH case 
while the blue dashed line is for the IH case. The horizontal, green long dashed lines show the 
region for the observed baryon asymmetry. Here we take M3 = 15 GeV, AMI2/MI = 10"*^, 
Recu = 7r/4, 6 = 77r/4, and r] = tt/3. We take sin^6'23 = 0.5 (left) and 0.36 (right), respectively. 
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A CP Asymmetry Parameter S^, 

In this appendix, we write down the CP asymmetry parameter 5^, defined in Eq. fl3ip . This 
parameter can be found from 



\^o.i\ a^ = ——^S,, (36) 

r=2..Sa=e.u.r \ / 



where we have neglected the terms which are proportional to AM/M^. 
In the NH case, we can parameterize S^, as 

6^ = di + rmd2 . (37) 

The leading term di is estimated as 

(ii = - sin 6^12 sin 2^^i3 [cos^ 6*13(3 + cos4^^23) — 4sin^ 613] sm{6 + rf) 

+ cos 6^12 sin 46^23 cos'^ 6'i3 sin 77 + 0{r^ , (38) 
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which had aheady shown in Eq. (132|) . For the sub- leading term we can find that 

(^2 = —2 sin 6^12 cos^ 6^12 sin 46*23 cos 6*13 (l + sin^ 6*12 cos 26*13) sin V 
— sin 6*12 cos^ 6*12 sin^ 6*23 sin 26*13 sin {6 — rj) 



+ i sin 2^13 



[-2 sin^ 6*12 cos 26*13 (7 + cos 46*23) + sin^ 26*23 (sin 6*12 + 5 sin 86*12)] sin {S + r]) 

o 

+2 sin^ 6*12 cos^ 6*12 sin 46*34 sin^ 6*13 cos 6*13 sin {26 + r]) . (39) 
In the IH case, on the other hand, we can parameterize S^, as 



6^ = (1 + r^)2 rfi + ^2, 



(40) 



where 

di = 



- sin 26*12 cos^ 



013 {-5 - 3 COS 4^23 + COS 2^13 (7 + COS 4^23)} 

— \ sin 4^12 {29 + 27 cos 4^23 + cos 4^13 (7 + cos 4^23) + 40 sin^ ^23 cos 2^13 1 
64 ^ ^ 

+ - sin 46*23 4 sin 6^13 {cos^ 6*13 + cos 6*12 (l + sin^ 6^13) } sin S cos?] 
8 L 

— {4 cos 26*12 sin 6*13 cos^ 6*13 + cos 46*i2 (7 sin 6*13 — sin 36*13) } cos 6 sin r] 

— - sin^ 2^23 sin^ ^13 (2 sin 26*i2 sin 26 cos 1] — sin 46*i2 cos 26 sin rj) , 



sin?7 



(41) 



d2 = 



^ sin 26*12 cos^ 6*13 {-5-3 cos 46*23 + cos 26*13 (7 + cos 46*23)} 
8 

sin 46*12 {29 + 27 cos 46*23 + cos46*i3 (7 + cos 46*23) + 40 sin^ 6*23 cos 26*13} 

64 



sin?7 



+^ sin 46*23 

o 



{sin 6^13 (1 — 7 cos 26^12) + 2 cos^ 6*12 sin 36*13} sin 6 cosr] 



— {4 cos 26*12 sin 6*13 cos^ ^13 + cos 4^12 



—7 sin 6*13 + sin 36*13) } cos 6 sin 77 



+ - sin^ 6*13 [sin 26*12 sin^ 26*23 sin 26 cos r] — sin 46*23 (1 — cos 46*23) cos 26 sin 77] . 
Actually, the leading term of Eq. f HOj) . i.e., 6i,\r^=o is obtained from di + c?2 as, 

= - sin 26*12 cos^ 6^13 {—5 — 3 cos 46*23 + cos 26*13 (7 + cos 46*23)} sin r] 



(42) 



di + d2 = - sin 26^12 

+4 sin 6*13 sin 46*23 (sin 6 cos r] 



{—5 — 3 cos 46*23 + cos 

— cos 26*12 cos 6 sin rj) 



This conclusion is consistent with Eq. 



(43) 
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